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Two-dimensional superconductors with time-reversal symmetry have a Z2 topological invariant, 
that distinguishes phases with and without helical Majorana edge states. We study the topological 
phase transition in a class-DIII network model, and show that it is associated with a metal-insulator 
transition for the thermal conductance of the helical superconductor. The localization length di- 
verges at the transition with critical exponent v » 2.0, about twice the known value in a chiral 
superconductor. 
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I. INTRODUCTION 

Gapped electronic systems are topological states of 
matter, subject to phase transitions in which a topo- 
logical invariant Q changes its valued For noninteract- 
ing electrons in two spatial dimensions, Q £ Z is inte- 
ger valued in the absence of time-reversal symmetry, and 
Q € I11 is binary in its presence. Two familiar examples 
from semiconductor physics are quantum Hall (QH) insu- 
lators in a strong magnetic field and quantum spin-Hall 
(QSH) insulators in zero magnetic fields The topologi- 
cal invariant Q determines the number of electrically con- 
ducting edge states, which changes by ±1 at a topological 
phase transition. The edge states arc chiral (unidirec- 
tional) in the QH effect and helical (counterpropagating) 
in the QSH effect. 

Superconductors can also have an excitation gap, with 
topologically protected edge states in two dimensions. 
The edge states carry heat but no charge, so supercon- 
ducting analogues of the quantum (spin)-Hall effect refer 
to the thermal rather than the electrical conductance. 
The thermal quantum Hall effect (2 6 Z) appears in 
the absence of time-reversal symmetry, for example, in a 
single layer of a chiral p-wave superconductor^— Time- 
reversal symmetry is restored in two layers of opposite 
chirality p x ± ip y , producing the thermal analogue of the 
quantum-spin Hall effect (Q £ Z2)£~— The edge states 
in both effects are Majorana fermions, chiral for the ther- 
mal QH effect and helical for the thermal QSH effect. 

There is a large literature on topological phase tran- 
sitions in QH and QSH insulators, as well as in chiral 
superconductors,— but the thermal QSH effect in helical 
superconductors has remained largely unexplored. This 
is symmetry class Dili, characterized by the absence of 
spin-rotation symmetry and the presence of both time- 
reversal and electron-hole symmetry. Here we present a 
study of the phase diagram and critical behavior in he- 
lical superconductors. We use a network model in sym- 
metry class Dili for a numerically efficient approach. 

We find that the main qualitative effect of time-reversal 
symmetry is that the transition between two topologi- 



cally distinct thermal insulators goes via a thermal metal 
phase, for finite but arbitrarily weak disorder. In con- 
trast, without time-reversal symmetry (in class D) the 
value of Q changes directly without an intermediate 
metallic phase for weak disorder. For strong disorder 
both chiral and helical superconductors have a thermal 
metal-insulator transition, but the critical behavior is dif- 
ferent: We find a localization length exponent v rj 2.0, 
about twice as large as the known value for chiral p-wave 
superconductors 

The outline of this paper is as follows. To put our re- 
sults for helical superconductors in the proper context, in 
the next section we first summarize known results for chi- 
ral p-wave superconductors. In Sec. IIIII wc introduce the 
network model of a helical superconductor, constructed 
out of two coupled chiral networks^ To identify topolog- 
ically distinct phases we apply a scattering formulation of 
the topological quantum number,— as described in Sec. 
IIVI In Sec. [V] we then present the results of our investiga- 
tion: the phase diagram with the thermal metal-insulator 
transition, the scaling of the thermal conductivity at the 
transition, and the critical exponent for the diverging lo- 
calization length. We conclude in Sec. EU 



II. CHIRAL VERSUS HELICAL TOPOLOGICAL 
SUPERCONDUCTORS 



According to the Altland-Zirnbaucr classification^ 
superconductors without spin-rotation symmetry are in 
class D or Dili depending on whether time-reversal sym- 
metry is broken or not. In two dimensions both symme- 
try classes can be in thermally insulating phases which 
are topologically distinct (with or without edge states). 
In this section we summarize what is known for the phase 
diagram in class D, before turning to the effects of time- 
reversal symmetry in class DHL 

A simple model Hamiltonian in class D represents a 
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FIG. 1: Solid curves (red): Phase boundaries of a chiral p- 
wave superconductor, symmetry class D, calculated in Ref. \13l 
from a lattice model based on the Hamiltonian (|2.ip .— The 
system is a thermal metal (Majorana metal) for strong dis- 
order and a thermal insulator for weak disorder. The dashed 
(green) lines show qualitatively the effect of time-reversal 
symmetry in class DHL the transition between two topologi- 
cally distinct insulators then goes via an intervening metallic 
phase. The critical exponents for the various phase transi- 
tions are indicated: the value for un is exact,— the values for 
vmi in class D (red, from Ref. Il3h and in class Dili (green, 
from this work) are numerical estimates. 



chiral p-wave superconductor in the x-y plane, 

H D = Va(PxT x +PyTy) + + U - fj)j t z . (2.1) 

The Pauli matrices tj (with tq the 2x2 unit matrix) 
operate on the electron-hole degree of freedom, coupled 
by the pair potential v/\(p x ±ip y ). The Fermi energy is /i 
and U(x, y) describes a random disorder potential (zero 
average) . 

By adding a spin degree of freedom (with Pauli ma- 
trices <Ji and unit matrix <jq), one can extend to the 
Hamiltonian of a helical superconductor in class Dili, 

Hmn = va(PxT x ct z + p v T y a ) + ( |~ + ^ ~ ^ ) T * a o 



KTyfJy. 



(2.2) 



Both Hamiltonians have electron-hole symmetry, 
t x H*t x = —H, but only the Hamiltonian (|2.2p has 
time-reversal symmetry, cfyH^ m a y = Hum- The term 
KTyUy in i?Din couples the two spin directions to zcroth 
order in momentum. Higher-order terms, such as pr z a y , 
can be included as well. 

The phase diagram of the chiral p-wave superconduc- 
tor was calculated for a lattice model in Ref. [lj, by dis- 
cretizing the Hamiltonian (|2.1j) on a lattice. A similar 
phase diagram was obtained earlier— ~— for a class-D net- 
work model (Cho-Fisher model).— As shown in Fig. [TJ 
there are two insulating phases plus a metallic phase at 
strong disorder. The two insulating phases are topolog- 
ical^ distinct, one is with and the other without chiral 



edge states. The disorder-induced thermal metal (Ma- 
jorana metal) arises because of resonant transmission 
through zero modes (Majorana fermions), pinned to po- 
tential fluctuations where U changes sign^Lr— 

The I-I phase boundary separating the two insulat- 
ing phases and M-I phase boundary separating insulat- 
ing and metallic phases meet at a tri-critical point. In 
the insulating phases the thermal conductivity decays 
exponentially cx e~ L ^ with system size L. The local- 
ization length £ diverges cx \(i — (i c \~ u on approach- 
ing a phase boundary at jj, = [i c . The critical expo- 
nent on the I-I phase boundary (at fi c = 0) is known 
analytically , 12 ' 24 v\\ = 1. The numerically obtained 
valued vya = 1.02 ± 0.06 on the M-I phase boundary 
is very close to v\\. Indeed, one would expect— z/mi = v \\ 
if the phase boundaries at the tricritical point meet at 
nonzero angle, which they seem to do. 

So much for a summary of known results for chiral su- 
perconductors in class D, without time-reversal symme- 
try. The time-reversally symmetric Hamiltonian -f/Din hi 
(|2.2p is just two uncoupled copies of H-q if K = 0. Upon 
increasing the coupling strength K, the time-reversal 
symmetry starts to qualitatively modify the phase dia- 
gram. As we will show in what follows (as is indicated 
schematically in Fig. [TJ, a metallic phase develops in 
between the two insulating phases at weak disorder for 
K £ 0. 



III. CLASS Dili NETWORK MODEL 

For numerical efficiency we use a network represen- 
tation of the class Dili Hamiltonian (|2.2p . Network 
models 2 ^ exist for quantum (spin)-Hall insulator o 27 ' 28 
and for chiral superconductors i 14 ' 29 A network model for 
helical superconductors was still missing and is provided 
here. 



A. Construction 

The network is defined on a two-dimensional bipartite 
square lattice, see Fig. [2] Helical Majorana modes propa- 
gate along the bonds and are scattered at the nodes. The 
hclicity of the modes signifies that the direction of motion 
is tied to the spin degree of freedom, f and J,, represented 
by dashed and solid lines in the figure. The modes en- 
circle local maxima and minima of the electrostatic po- 
tential (indicated by ± in the figure), in a clockwise or 
counterclockwise direction depending on the spin. The 
Majorana character of the modes signifies that there is 
no separate electron or hole mode, but one single charge- 
neutral mode per spin direction. 

The nodes of the lattice are saddle points between the 
local potential maxima and minima, alternating between 
adjacent plaquettcs in a checkerboard pattern. Scatter- 
ing at the nodes is described by 4 x 4 unitary scattering 
matrices S and S' that alternate between adjacent nodes 
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which amounts to a 90° rotation. 

The angle a characterizes the scattering at a node for 
each of the two spin directions separately, while the an- 
gle i? couples them. For i9 = we have two indepen- 
dent, time-reversed, copies of the Cho-Fisher modelp^ 
representing a pair of uncoupled chiral superconductors 
of opposite chirality. The orthogonal matrix O couples 
the two copies and produces a network model for a he- 
lical superconductor, in much the same way that Obusc 
et al.— obtained a network model for a quantum spin- 
Hall insulator by coupling a pair of Chalkcr-Coddington 
models^ for quantum Hall insulators in opposite mag- 
netic fields. 
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FIG. 2: Illustration of the network model described in the 
text. 

(black and white dots in the figure). The amplitudes 
o-na, b na of incoming and outgoing modes are related by 




= S 




= S' 




(3.1) 



(The labels a =f, I and n = 1, 2, 3,4 refer to Fig. |2J) 

A class Dili scattering matrix S at zero excitation en- 
ergy is constrained by both particle-hole symmetry and 
time-reversal symmetry In the basis (|3.ip (which relates 
time-reversed Majorana modes) the unitarity and sym- 
metry constraints read^ 



S = S* 



S , S — 1 . 



(3.2) 



so the scattering matrix is real orthogonal and antisym- 
metric. (The superscript T indicates the transpose.) 

The most general parameterization contains two real 
angles a, $ € (0, 27r) and one Z 2 index r) € { + 1, — l}r^ 



S 



o = 



A cos a 
O sin a 

- COS J? 

sin?9 



-O t sin a 
-7] A cos a 

-77 sin $ 
-77 cos $ 



A = 



1 

-1 



(3.3) 
(3.4) 



The corresponding parameterization for S' is obtained 
upon a permutation of the basis states, 



S' = PSP T , with P = 



'Q 1 0\ 
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(3.5) 



B. Vortices 

One difference between the superconducting network 
model considered here and the insulating model of Rcf. 
l28l . is that here the coupling of time-reversed networks 
is via real orthogonal rather than complex unitary ma- 
trices. This difference expresses the Majorana character 
of the modes, which have real rather than complex wave 
amplitudes. Another difference is the appearance of the 
I12 index 77, which determines the parity of the number 
of (time-reversally invariant) vortices in a plaquette. 

To see this, we take a — or a — ir, when the net- 
work consists of isolated plaquettcs. Denoting the value 
of i] for S and S' by rjs and r]s>, the phase factor ac- 
quired by the Majorana mode as it encircles a plaquette 
is r]sr]s>- A bound state at zero excitation energy (doubly 
degenerate because of time-reversal symmetry) results if 
VsVS' = 1- Since such Majorana zero- modes come in 
pairs, and each time-reversally invariant vortex in a he- 
lical superconductor traps one (doubly-degenerate) zero 
mode^ we conclude that the Majorana mode encircles 
an odd number of vortices for r)sr]s' = 1. 

In what follows we will assume that the system con- 
tains no vortices at all in the absence of disorder, so we 
choose r/s = 1, T)s' = — 1. The scattering matrices then 
take the form 



S = 



5" 





— r 

-i COST? 
t sin$ 

-tcostf 

tcostf 

r — £sin$ 
d, sini? 



icosi? — tsini?> 




(3.6a) 



(3.6b) 



-t COST? 



where we have abbreviated r = cos a, t = sin a. 



C. Vortex disorder 

Disorder is introduced in the network model by varying 
the scattering parameters in a random way from one node 
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to the other. We choose to keep the coupling strength 
•& the same for each node and to vary a. Following the 
same procedure as for the Cho-Fisher model^ we draw 
at at each node i independently from a distribution P{a.i) 
given by 

V{ai) = (1 - q)S(a l - a) + \q8{oii + a) 

+ ^qS(ai + a — tt). (3-7) 

The parameter q £ [0, 1] plays the role of disorder 
strength. 

This is a form of vortex disorder^ with probability 
q two time-reversally invariant vortices are inserted in 
the plaquettes adjacent to the i-th node, one vortex in 
one plaquette and another one in the diagonally opposite 
plaquctte. The diagonally opposite plaquettes are them- 
selves chosen with equal probability 1/2 from the two ± 
sublattices in Fig. [5] 

We use vortex disorder instead of purely electrostatic 
disorder [as in the Hamiltonians (|2.1[) and (|2.2[) ]. because 
it scatters more effectively and allows us to localize wave 
functions in smaller systems. Both forms of disorder can 
produce Majorana zero-modes ; 21 ' 33 so we do not expect 
qualitatively different features. 



IV. TOPOLOGICAL QUANTUM NUMBER 
AND THERMAL CONDUCTANCE 

A class Dili topological superconductor in two dimen- 
sions has a Z2 topological quantum number Q = ±1. 
Formulas for Q exist based on the Hamiltonian~i or on 
the scattering matrix^ Since the network model is de- 
scribed in terms of a scattering matrix, we use the latter 
formulation. 

We consider a rectangular geometry in the x-y plane, 
of length L = Af'a in the ^-direction and width W = Afa 
in the y-direction, where a is the lattice constant and 
Af, Af' are even integers (see Fig. [2]). In the trans- 
verse direction we impose either periodic boundary condi- 
tions, ip{x,Q) = ip(x,W), or antiperiodic boundary con- 
ditions, ip(x, 0) = —ip(x,W), on the wave functions. In 
the longitudinal direction we have absorbing boundary 
conditions^ corresponding to normal-metal reservoirs at 
x = and x = L. 

The scattering matrix S of the entire structure has 
Af x M reflection and transmission blocks, 

S = (* -£) ■ (4.1) 

The reflection block is a real antisymmetric matrix. Its 
Pfafhan 3 ^ determines the topological quantum number,— 

Q = signQ, Q= (Pf7e pbc )(Pfft apbc ), (4.2) 

where the labels pbc and apbc refer to the periodic and 
antiperiodic boundary conditions. (It does not matter if 
one takes 1Z or TV , they give the same Q.) 
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FIG. 3: Left panel: Phase diagram of the Dili network model 
without disorder. The solid curves, given by Eq. (|5.1[) . sep- 
arate helical topological superconductors with different val- 
ues Q of the Z2 topological quantum number. The red dot 
(at 1? = 0, 01 — 7r/4) marks the critical point in the Cho- 
Fisher mode^ of a chiral topological superconductor (class 
D). Right panel: Phase diagram in the presence of disorder 
(q = 0.07), when the topologically distinct phases are sepa- 
rated by a thermal (Majorana) metal. 

The use of periodic (or antiperiodic) boundary condi- 
tions is convenient to minimize finite-size effects and is 
sufficient to study bulk properties. To study edge prop- 
erties, one can impose reflecting boundary conditions by 
terminating the lattice as described in Ref. Ei?]. Depend- 
ing on how the lattice is terminated, one would then 
find that either Q = +1 or Q = - 1 produces a helical 
edge state along the boundary, and therefore represents 
a topologically nontrivial phase. In the present study, 
without reflecting boundaries, we can distinguish differ- 
ent topological phases — but we cannot decide which 
phase is trivial and which nontrivial. 

In addition to the topological quantum number we cal- 
culate the two-terminal thermal conductance G of the 
strip. This transport property is determined by the 
transmission matrix T, 

G = G TrTT t , (4.3) 

with Go = 7r 2 fc 2 jTo/6/i the thermal conductance quantum 
and To the temperature of the normal-metal reservoirs. 
The dimensionless thermal conductivity g is defined by 

g=(L/W)(G/G ). (4.4) 

For the calculation of the thermal conductivity we take 
periodic boundary conditions in the y-direction and a 
large aspect ratio W/L = 4, so that the choice of bound- 
ary conditions in the transverse direction has only a small 
effect. 



V. TOPOLOGICAL PHASE TRANSITIONS 

A. Phase diagram without disorder 

In the absence of disorder the location of the topolog- 
ical phase transitions can be determined exactly using 
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FIG. 4: Data points (with dashed line as a guide to the eye): 
conductivity as a function of system size (at fixed aspect ratio 
W/L = 4) of the Dili network model without disorder, at 
criticality for ■& = 0.4. The ballistic scaling g oc L results from 
the Fermi circle of the long- wavelength Hamiltonian ()5.2|) (red 
dotted circle in the inset). 



a duality relation, see App. [A] We find that the topo- 
logical quantum number switches sign at critical points 
a c ,'d c that satisfy 



sin a r cos d r 



(5.1) 



The phase boundaries are plotted in Fig. [3] (left panel) , 
together with the values of the topological quantum num- 
ber (|4.2p at the two sides of the transition. For d c = we 
recover the known value a c = 7r/4 of the critical point 
in the chiral Cho-Fishcr models This is as expected, 
since for & = our helical network model consists of two 
independent chiral copies. 



B. Scaling of the critical conductivity 



At the phase boundaries the excitation gap of the sys- 
tem closes. In the chiral network model this produces 
a scale-invariant thermal conductivity g = l/vr, regard- 
less of whether the system is disordered or not^ but in 
our helical model the conductivity at the critical point 
scales ballistically oc L in the absence of disorder. The 
ballistic scaling is demonstrated in Fig. |4] and can be un- 
derstood by examining the long- wave length Hamiltonian 
corresponding to the network model. 

The calculation of this Hamiltonian proceeds entirely 
along the lines of the analogous calculation for the quan- 
tum spin-Hall insulator^ so we only give the result. To 
first order in the deviation from the Cho-Fishcr critical 
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FIG. 5: Disorder average of Q, defined in Eq. (|4.2|) as a func- 
tion of q for fixed d = 0.5, q — 0.1, and different values of 
L = W/4. The topological quantum number switches from 
— 1 to +1 in the insulating phases via a plateau of zero average 
in the metallic phase. 



point (a,i?) = (7r/4, 0), we find 
H = [p x r x a z + PyTyO-Q + 2(a - 7r/4) r z a 

= {PxT x +PyT y )a Q + 2(a - 7r/4) t z <j z + V2-dr a x 



U 



t q + ir y )a 



(r - ir y )a z . 



(5.2) 
(5.3) 



Up to a unitary transformation U, and to first order in 
(dimcnsionless) momentum p, this Hamiltonian has the 
form of Hdui in Eq. (|2.2I) . with a playing the role of the 
chemical potential [i and $ playing the role of the spin 
coupling strength K . 

The gap closes at a = 7t/4 on a (twofold degenerate) 
circle p x +Py = 2i9 2 in momentum space (inset in Fig. 
2]). The Fermi wavevector kp = \/2|i?|/a corresponds to 
a ballistic conductance G/Gq = 2kpW ~/ir. Hence we find 
the ballistic critical conductivity 



9c 



2y/2\i}\L 



< 1. 



(5.4) 



The ballistic scaling of the critical conductivity is a sig- 
nature of the appearance of a Fermi circle at the phase 
transition, which is a special property of our model (cho- 
sen to maximize the coupling between the two oppo- 
site chiralitics) . More generally, the gap in a class-DIII 
Hamiltonian will close at four isolated points in momen- 
tum spaced resulting in a scale-invariant critical con- 
ductivity. 



C. Phase diagram with disorder 

As disorder is introduced in the system, a metallic 
phase develops in between the insulating phases, so that 
the switch from Q = +1 to Q = — 1 goes via two metal- 
insulator transitions. In the metallic region Q has a ran- 
dom sign, averaging out to zero (see Fig. [51 where we 
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FIG. 6: Phase diagram of the disordered Dili network model 
for fixed ■& = 0.785. 




FIG. 7: Scaling of the conductivity in the disordered Dili 
network model (a = 1.2, $ = 0.5), at the metal-insulator 
transition (blue data points, q — 0.175) and in the metallic 
phase (red data points, q = 0.2). The dotted lines indicate 
the scale-invariance of the critical conductivity and the loga- 
rithmic scaling of the metallic conductivity. 



averaged Q rather than Q = sign Q to reduce statisti- 
cal fluctuations) . The phase diagram is shown in Fig. [3] 
(right panel) for a fixed disorder strength q in the a-d 
plane and in Fig. [5] for fixed $ in the a-q plane. The 
metallic regions become broader and broader with in- 
creasing disorder, and for q > 0.2 no insulating phase is 
left. 

Fig. [7] shows the sample-size dependence of the av- 
erage conductivity, both at the metal-insulator transi- 
tion and in the metallic phase. (The exponential de- 
cay oc exp(— L/£) in the insulating phase is not shown.) 
While without disorder the conductivity scales ballisti- 
cally oc I at the critical point (see Fig. 2]), disorder re- 
stores the scale invariance that is the hallmark of criti- 
cality. In the metallic phase we find a logarithmically in- 
creasing conductivity (g) = cln(L/a), characteristic for a 
Majorana metal^ ' 12 i 13 , with c = 1/tt (dotted line in Fig. 
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FIG. 8: Disorder-averaged conductivity as a function of dis- 
order strength for various system sizes L = W/A at fixed 
a — 1.2, ■& = 0.5. Solid curves are a fit to the scaling law, as 
described in App.[B] The curves cross at the critical point. 



D. Critical exponent 

The metal-insulator transition is associated with a di- 
verging localization length £ oc \x— x c \~ v , where x can be 
any of the control parameters a, q and x c is the value 
of x at the critical point. To determine the critical ex- 
ponent v we perform a finite-size scaling analysis of the 
thermal conductivity, in a manner analogous to the work 
of Slevin and Ohtsuki4£ 

Typical data is shown in Fig. [8j where we follow the 
thermal conductivity through the metal-insulator transi- 
tion upon varying the disorder strength q at fixed a and 
■d. The dimensions L and W of the system are increased 
at fixed aspect ratio W/L = 4. The curves are fits of the 
data to the scaling law, as described in App. [B] Results 
are given in Table |TJ 



control parameter 


V 


fjc 


x = q, a = 1.2, i9 = 0.5 


2.06 [1.89,2.20] 


0.74 [0.66,0.81] 


x ee ct, q = 0.1, = 0.785 


1.93 [1.78,2.24] 


0.73 [0.69, 0.76] 



ED. 



TABLE I: Critical exponent v and critical conductivity g c , 
with ±3cr confidence levels indicated by square brackets, ex- 
tracted from a finite-size scaling analysis with x = q or x — a 
as control parameter to drive the system through a metal- 
insulator transition. The two cases correspond to transitions 
into a topological phase with opposite value of Q. 



VI. CONCLUSION 

In conclusion, we have presented a network model for a 
two-dimensional helical p-wave superconductor and used 
it to investigate the topological phase transitions. The 
presence of time-reversal symmetry (class Dili) leads to 
differences with the more familiar chiral p-wave super- 
conductors (class D), of a qualitative nature (the ap- 
pearance of a thermal metal separating the thermally 
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TABLE II: Overview of the critical exponents in the five symmetry classes that exhibit a topological phase transition in two 
dimensions. 



insulating phases) and also of a quantitative nature (an 
approximate doubling of the critical exponent). Helical 
superconductors have not yet been convincingly demon- 
strated in experiments, but there is an active search and 
a variety of candidate materials 42r— 

This study fills in the last missing entry in the list of 
critical exponents of two-dimensional topological phase 
transitions (see Table HI)) , completing a line of research 
on network models that started with the seminal work of 
Chalker and Coddington on the quantum Hall effect^ It 
is intriguing that the effect of time-reversal symmetry is 
close to a doubling of the critical exponent (from v w 1.0 
in class D to v « 2.0 in class Dili), but this may well be 
accidental. 
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Appendix A: Location of the critical point in the 
network model without disorder 

The critical point in the clean Dili network model 
(q = 0) can be obtained from a duality relation: Ex- 
change of the scattering matrices S and 5" of the two 
sublattices (black and white nodes in Fig. [5]) has the ef- 
fect of exchanging the trivial and nontrivial phases. This 
can be seen most easily for reflecting boundary condi- 
tions, when exchange of the sublattices either creates or 
removes the helical edge state (see Fig. 7 of Ref. H3) ■ 

The exchange of S and S' amounts to the transforma- 
tion of a, $ into a', given by 

cos a' = — sin a cos , cos a = — sin a' cos , 

sin a' sin $ = sin a sin i?. (Al) 

Equivalently, the unit vector h = 
(sin a cos sin a sin cos a) is transformed into 

n' x = -n z , n' y = n y , n' z = -n x , (A2) 



which amounts to a reflection in the plane x + z = 0. The 
network is selfdual if n lies in this plane. Since a selfdual 
network is at the critical point, we arrive at a sufficient 
condition for criticality, n x + n z =0, or equivalently 

sina c cosi9 c + cosa c = 0. (A3) 

An alternative condition can be obtained by noting 
that the transformation a t— > — a leaves the reflection 
matrix 1Z unaffected. The topological quantum number 
(|4.2p therefore remains unchanged, so if a Cl -d c is criti- 
cal then also —a Cl d c . We thus have a second sufficient 
condition for criticality, 

— sina c cos^ c + cosa c = 0. (A4) 

Eqs. (|A"3|l and (|A~4|) together give the condition (|5T|) . 



Appendix B: Finite-size scaling analysis 

We determine the critical exponent v associated with 
the metal- insulator transitions on both sides of the metal- 
lic phase, by an analysis of the system size L = W/4 
dependence of the disorder averaged conductivity g. Fol- 
lowing the general approach of Slevin and Ohtsuki^S we 
take into account finite-size corrections to scaling in the 
form of nonlincaritics in the scaling variable u, as well as 
the presence of an irrelevant scaling exponent y < 0. 

The finite-size scaling law reads 

g = F(u L l '\u 1 Ly) 1 (Bl) 

in terms of the relevant scaling variable uq and the lead- 
ing irrelevant scaling variable u\. We perform a Taylor 
expansion, first up to order n in powers of u\ , 

n 

g = Y / 4L ky Fk(u L 1 /»), (B2) 

fe=0 

and then on each of the functions Fk up to order rrik in 
powers of uq, 

F k {u Q L 1 ' v ) = Y,niu' v F ki . (B3) 
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control 
parameter 


fit parameters 


irrelevant exp. 


reduced 

x 2 


nr. of degrees 
of freedom 


goodness 
of fit 


x = q 
a = 1.2, i? = 0.5 


n = mi = 1, mo = 3 
q r = 2, qi = 


y = -0.44 
[-0.81,-0.30] 


1.04 


97 


0.65 


x = a 
g = 0.1, 1? = 0.785 


n = 1, mo = mi = 2 
g r = 3, g» = 1 


y = -0.67 
[-0.78, -0.53] 


1.06 


128 


0.56 



TABLE III: Parameters for the nonlinear fitting analysis, giving the critical exponent and conductivity of Table HT1 



We tune through the metal-insulator transition by 
varying one parameter x £ {a, q} through the critical 
point x c , keeping the other two parameters fixed. Non- 
linearities are taken into account by Taylor expanding 
the relevant and irrelevant scaling variables in powers of 
x — x c , up to orders q r and qi respectively, 

q r 

u (x-x c ) = ^b k (x-x c ) k , (B4) 

k=l 

ui{x - x c ) = c k (x - x c ) k . (B5) 

fc=0 



The expansion of the relevant scaling variable does not 
contain a zcroth order term, due to the requirement 
uo(0) =0 for a scale-invariant critical conductivity. 

The average conductivity is determined up to a pre- 
cision between ~ 0.2% and ~ 0.07% (error bars much 
smaller than the size of the symbols in Fig. We per- 
form the fit by minimizing the % 2 statistic, and express 
the goodness of fit as well as the degree of uncertainty 
in the fit parameters through a Monte Carlo resampling 
technique^ as appropriate for a non-linear fitting func- 
tion. Results are collected in Tables Q] and IIIII 
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